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i—^ , Abstract. Modelling amplitudes of stochastically excited oscillations in stars is a powerful tool 

C/j ' for understanding the properties of the convective zones. For instance, it gives us information 

on the way turbulent eddies are temporally correlated in a very large Reynolds number regime. 
We discuss the way the time correlation between eddies is modelled and we present recent 
theoretical developments as well as observational results. Eventually, we discuss the physical 
^ \ underlying meaning of the results by introducing the Ornstein-Uhlenbeck process, which is a 

■ sub-class of a Gaussian Markov process. 
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Amplitudes of solar-like oscillations result from a balance between excitation and damping and 
crucially depend on the way the eddies are temporally correlated (see the review of [1]). Hence, 
the improvement of our understanding and modeling of the temporal correlation of turbulent 
eddies, hereafter denoted in the Fourier domain as Xki^), is fundamental to infer turbulent 
properties in stellar convection zones. 

Most of the theoretical formulations of mode excitation explicitly or implicitly assume a 
Gaussian functional form for Xki^) ([2^^]). However, 3D hydrodynamical simulations of the 
outer layers of the Sun show that, at the length-scales close to that of the energy bearing eddies 
(around 1 Mm), Xk is a Lorentzian function [8, 9]. As pointed-out by [7], a Lorentzian Xk is also 
a result predicted for the largest, most-energetic eddies by the time-dependent mixing- length 
formulation derived by [10]. However, [7], [1], and [11] found that a Lorentzian Xk^ when used 
with a mixing-length description of the whole convection zone, results in a severe over-estimation 
for the low-frequency modes. 

In a recent work, [12] introduced a cut-off frequency in the modelling of the eddy-time 
correlation function to account for the effect of short-time scales. Indeed, under the sweeping 
approximation, which consists in assuming that the temporal correlation of the eddies, (in the 
inertial subrange) is dominated by the advection by energy-bearing eddies, the shape of the 
temporal correlation function of eddies is no longer Lorentzian. 

Hence, in this paper we discuss the modelling of the eddy-time correlation function in both 
the limit of large and short time scales. This modelling is then validated using observational 
constraints from amplitudes of solar-like oscillations. Eventually, the underlying physical 



meaning of the observed shape of the temporal correlation function is discussed by mean of 
a simplified stochastic process. 



2. Modelling the Eulerian eddy time-correlation function 

The formalism we used to compute excitation rates of radial modes was developed by [13] (see 
also [14]). For a radial mode of frequency ujq = 2ttuo, the excitation rate (or equivalently, the 
energy injection rate), P, mostly arises from the Reynolds stresses and can be written as (see 
Eq. (21) of [15]) 
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where is the radial component of the fluid displacement eigenfunction (S^), m is the local 
mass, pq the mean density, ujq the mode angular frequency, / the mode inertia, Sk the source 
function, E{k) the spatial kinetic energy spectrum, Xk the eddy-time correlation function, and 
k the wave- number. 



2.1. Some preliminary definitions 

We now turn to a rigourous definition of the eddy-time correlation function (xfc)- For a turbulent 
fluid, one defines the Eulerian eddy time-correlation function as 

\k x 



{u{x + f,t + T)-u{x,t)) = y f (/k, t, r) e"'-^ d^/fc , (3) 

where u is the Eulerian turbulent velocity field, x and t the space and time position of the fluid 
element, k the wave number vector, r the time-correlation length, and r the space-correlation 
length. The function £ in the RHS of Eq. (3) represents the time-correlation function associated 
with an eddy of wave-number k. 

We assume an isotropic and stationary turbulence, accordingly £ is only a function of k and 
T. The quantity £{k,T) is related to the turbulent energy spectrum according to 

SiKr,^^- (4) 

where E{k, r) is the turbulent kinetic energy spectrum whose temporal Fourier transform is 
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E{k,uj) = — E{k,T)e''^^ dr (5) 

where oj is the eddy frequency, and E{k,uj) is written as follows [13, 16] 

r+oo 

E{k,u) = E{k)xk{^) with / Xk{^)du = l (6) 



where Xk{'^) is the frequency component of E{k,oj). In other words, Xki^) represents - in the 
frequency domain - the temporal correlation between eddies of wave-number k. 



2.2. Long-time scales 



A modelling of Xk in the stellar context from the full non-linear hydrodynamic equations is 
a very difficult task, and no such a model has been proposed so far. This is particularly the 
case for large times, which we defines by the times of the order or greater than the integral 
time-scale that is the width at half maximum of Xk- Hence, we adopt an alternative approach 
that is to postulate the functional form of Xk and then to confront it to observations and/or 3D 
hydrodynamical numerical simulations. 

Several functional have then been tested in the past and as mentioned in Sect. 1, theoretical 
and observational evidence show that Xk{^) is Lorentzian, i.e. 

Xk{^) = : — ; — -o with ujk = kuk and ul = E{k)dk (7) 

vrcjfc 1 + {oj/uJk) Jk 

where ujk is by definition the width at half maximum of Xki^^)- In the framework of [13] 's 
formalism, this latter quantity is evaluated as: 

ujk = kuk with 1^^= I E{k)dk (8) 

Jk 

where E{k) is defined by Eq. (6). Note that Eq. (7) corresponds to an exponential in the time 
domain. 

Such a modelling of the temporal correlation of eddies has successfully been tested in the solar 
and stellar context (see [1] for a detailed discussion). However, when coupled with a standard 
ID model of the solar convection zone, it leads to an important overestimation of the solar mode 
amplitudes. More precisely, a substantial fraction of the energy injected to those modes comes 
from the deep layers of the solar convective region. As a result, [7] and [1] then suggested that 
most contributing eddies situated deep in the Sun have a Xk rather Gaussian than Lorentzian 
since, at fixed frequency, a Gaussian Xk decreases more rapidly with depth. Nevertheless, this 
conclusion has recently been questioning by [12], which demonstrate that Eq. (7) is no longer 
valid for short-time scales. 

2.3. Short-time scales 

To infer some properties of the eddy-time correlation function at short-time scales, the function 
£{k, t, t) appearing in Eq. (3) is expanded for short time scales. In the inertial sub-range it gives 
(see [17] for a detailed derivation) 

£{k, t) = £{k, r = 0) (^1 - Qfclrl - i {ujErf + . . . ^ (9) 

where the characteristic frequency can be estimated by the eddy turn-over frequency cok 
(see [12] for details), and the second characteristic frequency, ujE{k), is the curvature of the 
correlation function near the origin, and defined by 

oje = kuo . (10) 

This expression is obtained under the sweeping assumption, which consists in assuming that 
the velocity field il{k) associated with an eddy of wave-number k lying in the inertial-subrange 
(i.e. large k compared to ko) is advected by the energy-bearing eddies with velocity uq {i.e. of 
wave- number ko). 
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Figure 1. Left panel: Solar p-mode excitation rates as a function of the frequency v. The 
dots correspond to the observational data obtained by the GONG network, as derived by [18], 
and the triangles corresponds to observational data obtained by the GONG network as derived 
by [19] for £ = to ^ = 35. The dashed line corresponds to the computation of the excitation 
rates using with a Lorentzien description of Xk without any cut-off frequency, while the dashed 
triple dot line corresponds to the same computation except that the Lorentzian description of 
Xk down to the cut-off frequency uje- The solid line corresponds to the computation of mode 
excitation rates using 3D numerical computation to described turbulent convection properties 
and a Lorentzian Xk together with a cut-off frequency at uj = uje- Note that no scaling is applied 
to those calculations. Right panel: Ratio between (Vmax) the maximum of the mode velocity 
relative to the observed solar value {VQ^max = 25.2cm/ s for £ = 1 modes, see Kjeldsen et al. 
2008). Filled dots correspond to the stars for which solar-like oscillations have been detected in 
Doppler velocity (see [20]. The lines correspond to the power laws obtained from the predicted 
scaling laws for Pmax and estimated values of the damping rates rjmax (see [1] for details). 



To go further, we note that the zeroth- and first-order terms in Eq. (9) are consistent with an 
exponential decrease of width Ofc (i.e. a Lorentzian in the frequency domain of width , Eq. (7)) 
for small r. In contrast, the zeroth-order term together with the second order term in Eq. (9) 
are consistent with a Gaussian behavior of width te- In turn, [12] have shown that for ui > loe 
the second order term dominates over the first order one in Eq. (9), at all length-scales. Hence, 
for frequencies near the micro-scale frequency (a; > uje), the eddy-time correlation function 
behaves as a Gaussian function (e~('^/'^^) ) rather than a Lorentzian function, resulting in a 
sharp decrease with uj. Hence, the contributions for uj > uje are negligible. 

3. Observational evidences for a Lorentzian eddy-time correlation function 
3.1. The solar p -mode energy injection rates 

When the frequency range of Xk is extended toward infinity, computation of P fails to reproduce 
the observations, in particular the low-frequency shape (see Fig. 1). It is in agreement with the 
results of [7] and [1] , and results in an over-estimation of the excitation rates at low frequency. 
In contrast, by assuming that the time-dynamic of eddies in the Eulerian point of view is 
dominated by the sweeping, the Eulerian time micro-scale arises as a cut-off frequency (see 
Sect. 2.3). Hence, Xki^) is modeled following Eq. (7) for uj < uje and Xki^) = elsewhere. 
Using such a procedure to model Xk{^) permits us to reproduce the low-frequency < 3mHz) 
shape of the mode excitation rates as observed by the GONG network (see Fig. 1). This is 



explained as follows; for large scale eddies near , situated deep in the convective region, the 
cut-off frequency w^; is close to uj^- As a consequence, the frequency range over which Xk is 
integrated in Eq. (2) is limited, resulting in lower injection rates into the modes. This results 
then re-conciliates theoretical and observational evidences that the frequency dependence of 
the eddy-time correlation must be Lorentzian in the whole solar convective region down to the 
cut-off frequency uje- 

Note that a modelling of the excitation rates, using a Gaussian for Xk and a 3D numerical 
simulation to model convection, leads to an underestimation of the excitation rates, as shown 



3.2. Scaling law on the maximum mode energy injection rate for solar-type pulsators 

A similar conclusion has been reached by [21] for main-sequence and red giant stars. Indeed, 
[22] have shown that the maximum mode surface velocity scales as a function of the ratio of 
luminosity to the mass {L/M). Later, [21] have shown that the exponent (3 is sensitve to the 
adopted model for computing Xfc; ^-e. /3 = 0.7 for a Lorentzian and (3 = 1 for a Gaussian. 
As shown in Fig. 1, the best agreement with the observations is found when a Lorentzian Xk is 
assumed. In contrast, assuming a Gaussian Xk results in a larger exponent (3 and do not permit to 
reproduce the observations. Consequently, we conclude that the eddy-time correlation function 
is a Lorenztian functional (or equivalently an exponential functional in the time domain) for all 
solar-like pulsators and in the upper-layers of convective regions. Note that the case of deep 
layers in convective regions is still to be investigated, but 3D numerical results from [9] also 
suggest Xk is Lorentzian. 

4. On the interpretation of the Lorentzian shape of Xk 

A modelling of the Eulerian time correlation function from the full non-linear hydrodynamic 
equations is a difficult task, making the interpretation of our results non-trivial. However, one 
can use a classical approach in turbulence, which consists in modelling the turbulence by a 
stochastic model (see the review [23]). Similarly to the Brownian motion, one can define a fluid 
elements and express the corresponding Langevin equation. Naturally, a Lagrangian approach 
is more adapted however, our modelling of mode amplitudes deals with Eulerian velocities. To 
overcome this difficulty we first introduced a relation between Eulerian and Lagrangian time 
correlation function, first introduced by [24]. 

^.1. The Corrsin's conjecture 

The problem of relating Eulerian and Lagrangian correlation is an outstanding problem in 
turbulence, since the measurement of Lagrangian quantities is experimentally difficult. A general 
expression can be written such as 



where w is the Lagrangian velocity and u the Eulerian velocity, and ip is the probability 
distribution function of the fluid particule displacement at time t. In this framework, [24] 
proposed a simplified relation based on the assumption that for large times (r), ip and u are 
statistically weakly dependent. Then by using a Gaussian distribution for ■0, [24] demonstrates 
that one can rewritte Eq. (11) in the limit of large scales^ 



by [12]. 







^ This assumption is valid for our purposes since modes are mainly excited by large scales. 



This approximated relation has been extensively discussed and has been shown to be rather 
accurate for turbulent flows without helicity (e.g. [25, 26]). 



4-2. The Ornstein-Uhlenbeck Process 

For long time-scales, given the relation Eq. (12), we now turn to discuss some aspect of 
Lagrangian stochastic models. In the limit of a stationary, homogeneous turbulence, a stochastic 
model can be adopted that follows the Langevin equation (written in term of stochastic 
differential equation) 

dw{t) = -w{t)dt/T + {2u'^ /T)Cit)dt (13) 



where T is the Lagrangian integral time-scale, n'^ is the variance of wit), and C a random 
fluctuating term. To go further, lets us assume that is a Gaussian random variable with zero 
mean and which is rapidly varying, i.e. < C(OC(^') >= ^{t ~ i')- 

Consequently, the Langevin equation Eq. (13) describes a continuous Gaussian Markov 
process in time (it is a Wiener process or equivalently a Brownian process). More precisely, 
since we are interested to describe a stationary process, w{t) is a stationary, Gaussian, Markov 
process know as the Ornstein-Uhlenbeck Process (see [27] for a more detailed definition). In 
other terms it corresponds to a noisy relaxation process. A general property of such a process 
is that its time correlation function is (e.g., [23, 27]) 
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which, in the time Fourier space, corresponds to a Lorentzian function. 

A more rigourous investigation has been proposed by [28], who proposed a stochastic model 
for turbulent dissipation. The author then proposed an expression for the time correlation 
function function of the Reynolds number 



{w{t + T)w{t)) 
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(15) 



which tends to Eq. (14) for Re — )• -|-oo. 

Hence, by used of the Corrsin's conjecture one can conclude that our results, obtained using 
mode amplitudes {i.e. Xk is Lorentzian in the time Fourier domain), suggests that turbulent 
eddies in convective region of solar-type stars follow such a noisy relaxation process. 



5. Conclusion 

We have shown that the modelling of mode amplitudes gives us access to some properties of 
turbulent convection in stars. More precisely, we detailed the way the time correlation between 
eddies is actually modelled and validated using observational data. The conclusion is that the 
eddy-time correlation function follows a Lorentzian shape (in the time Fourier domain) down to 
a cut-off frequency, in convective region of solar-type stars. 

We then discussed this shape and we emphasised that a Lorentzian function - or equivalently 
an exponential function in the time domain - is a rather common result for noisy relaxation 
processes. Our discussion also shows that maybe the more interesting will be to investigate the 
departure from a Lorentzian shape of Xk^ and its physical interpretation. To this end, we stress 
that more accurate observations for the Sun of solar-like pulsators is mandatory. 
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